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ABSTRACT 

In an earlier paper we quantified the mean merger rate of dark matter haloes as a 
function of redshift z, descendant halo mass Mg, and progenitor halo mass ratio ^ us- 
ing the Millennium simulation of the ACDM cosmology. Here we broaden that study 
and investigate the dependence of the merger rate of haloes on their surrounding envi- 
ronment. A number of local mass overdensity variables, both including and excluding 
the halo mass itself, are tested as measures of a halo's environment. The simple func- 
tional dependence on z, Mq, and ^ of the merger rate found in our earlier work is 
largely preserved in different environments, but we find that the overall amplitude of 
the merger rate has a strong positive correlation with the environmental densities. 
For galaxy-mass haloes, we find mergers to occur ~ 2.5 times more frequently in the 
densest regions than in voids at both z — and higher redshifts. Higher-mass haloes 
show similar trends. We present a fitting form for this environmental dependence that 
is a function of both mass and local density and is valid out to z = 2. The amplitude 
of the progenitor (or conditional) mass function shows a similar correlation with local 
overdensity, suggesting that the extended Press-Schechter model for halo growth needs 
to be modified to incorporate environmental effects. 



1 INTRODUCTION 



& Cole 1993 1 that are widely used for making theoretical 



In studies of cosmological structure formation, the mass of 
a dark matter halo is a key variable upon which many prop- 
erties of galaxies and their host haloes depend. For instance, 
dark matter haloes of lower mass are expected to form earlier 
on average than more massive haloes in hierarchical cosmo- 
logical models such as ACDM. In semi-analytical modelling 
of galaxy formation (see |Baugh||2006] for a review), proper- 
ties such as the formation redshift, halo occupation number, 
galaxy colour and morphology, and stellar vs AGN feedback 
processes are all assumed to depend on the mass of the halo 
(sometimes better characterised by the halo circular veloc- 

ity). 

In addition to the halo mass, however, recent work 
based on numerical simulations has shown that a halo's lo- 
cal environment also affects various aspects of halo forma- 
tion. For instance, at a fixed mass, older haloes are found 
to cluster more strongly than more recently formed haloes 
([Gottlo ber et al.|[200l| [Sheth fc Tormen| [2004l | Cao et al 
2005; Harker et al.| |2006| [Wechsler et aU ^OOe Ji ng et al 
2007| |Gao fc White] |2007, Maulbetsch 



predictions of galaxy statistics and for Monte Carlo con- 
structions of merger trees. The lack of environmental cor- 
relation arises from the Markovian nature of the random 
walks in the excursion set model. This limitation is not built 
into the model per se, but is an assumption stemming from 
the use of a tophat Fourier-space window function. When 
a Gaussian window function is used, for instance, IZentner] 
( |2007[ ) finds an environmental dependence in the halo for- 
mation redshift, but the dependence is opposite to that seen 
in the numerical simulations cited above. Other attempts at 
incorporating environmental effects into the excursion set 
model thus far have not been able to reproduce the corre- 



2007, ,Wang et al. 



et al. 20071. Other halo properties such as concentration, 
spin, shape, and substructure mass fraction have also been 



found to vary with halo environment (e.g., Avila- Reese et al. 
2005; Wechsler ot al. 20061 1 Jing et al.||2007| |Gao fc Whit~ 
'2007(|Bett et al.||2007| ). 



lations in simulations (e.g., Sandvik et al.|2007 Desjacques 
|2008j). 

In this paper, we focus on the environmental depen- 
dence of the merger rate of dark matter haloes, a topic that 
has not been studied in detail. The merger rate is an im- 
portant quantity for understanding and interpreting obser- 
vational data on galaxy formation, growth, and feedback 
processes. While the mergers of galaxies and the mergers 
of dark matter haloes are not identical processes, the two 
processes are closely related, and quantifying the latter is 
the first key step in understanding the former. There have 
been few theoretical studies of merger rates (e.g.,|Gottlober| 



In contrast, no such environmental dependence is pre- 
dicted in the extended Prcss-Schcchtor (EPS) and excursion 
set models ( Press fc Schechter|1974||Bond et al.|1991 Lacey 



et al.|2001||Fakhouri fc Ma|2008[ [Stewart et al.|2008[ ) prob^ 
ably because mergers are two- (or more-) body processes, 
and a large ensemble of descendent haloes and their pro- 
genitor haloes must be identified from merger trees before 
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the rate can be reliably calculated. In comparison, studies 
of halo properties such as the mass function, density and 
velocity profiles, concentration, triaxiality, spin, and sub- 
structure distribution require only the particle information 
from a single simulation output. 

This paper is an extension of our earlier study |Fakhouri| 
fc Ma| (|2008!) (henceforth FM08). There we quantified the 



global mean merger rates of haloes in the Millennium simu- 
lation ( Springel et al.|2005 1 over a wide range of descendant 



halo mass (10 < Mo < 10 M0), progenitor mass ratio 
(10"^ ^ C < 1), and redshift (0 z < 6). We found that 
when expressed in units of the mean number of mergers per 
halo per unit redshift, the merger rate has a very simple de- 
pendence on Mo, 5, and z: the rate depends very weakly on 
halo mass (oc A/q"®) and redshift, and scales as a power law 
in the progenitor mass ratio (oc ''^) for minor mergers 
(^ < 0.1), with a mild upturn for major mergers. These sim- 
ple trends allowed us to propose a universal fitting form for 
the mean merger rate that is accurate to 10-20%. 

Here we go beyond the global merger rate and use the 
rich halo statistics in the Millennium database to quantify 
the merger rate as a function of halo environment, in ad- 
dition to descendant mass, progenitor mass ratio, and red- 
shift. We also investigate the environmental dependence of 
the progenitor (or conditional) mass function. This quantity 
is closely related to the merger rate and is also the most 
important ingredient in the EPS and excursion set models 
for constructing Monte Carlo merger trees. 

Several recent environmental studies have used halo 
clustering, quantified by the halo bias, as a measure of envi- 



ronment (e.g. 


Gottlober et al.|2002| 


Sheth & Tormen||2004 


Gao et al.|2005 Harker et al. 2006||Jing et al.|2007||Wechsler| 


et al. 2006 


Gao & White 2007 


1 . While halo bias is a powerful 



statistical quantity, we choose a simpler and more intuitive 
local environment measure and use the local mass density 
centred at each halo. The earlier studies that have used local 
overdensities as measures of halo environment have used a 
variety of definitions, e.g., the mean density within a sphere 
of some radius (ranging from 4 to lQh~^ Mpc) o r within 
a spherical shell (e.g. between 2 and 5/i~^ Mpc) (Lemson 
fc Kauffmann|[T999l [Haxker et al.||2006l IWang et al.||2007| 



Maulbetsch et al.||2007| |Hahn et ai.||2008| . In this paper we 
compare different definitions of the local overdensity, both 
including and excluding the mass of the central halo itself. 

In this paper we also provide an in-depth investigation 
of the effects of halo fragmentation on the merger rate and 
its environmental dependence. In FM08, we discussed how 
fragmentation is a generic feature of all merger trees and 
compared the stitching method with the conventional snip- 
ping method for handling these events. We will show here 
that fragmentation occurs more frequently in dense regions 
than in voids; understanding the effects of fragmentation on 
merger rates is therefore essential for obtaining robust re- 
sults in dense environments. There are three general types 
of approaches to handling fragmentations: do nothing [snip- 
ping), stitching together fragmented haloes, or splitting up 
the common progenitor of the fragmented haloes. We will 
compare five algorithms for handling fragmentations based 
on these three approaches and show that except for one algo- 
rithm, all the algorithms give similar merger rates to within 
20%. 

This paper is organised as follows. In § [2] we briefly re- 



view how haloes and merger trees are constructed from the 
particle data in the Millennium simulation. Statistics detail- 
ing the distribution of halo mass at different redshifts are 
summarised in Table [l] In §|3]we compare four local density 
measures and their distributions in relation to halo mass. 
Three of the measures use the dark matter mass in a sphere 
centred at a given halo, either including or excluding the 
mass of the central halo. The fourth measure is motivated 
by observables such as luminosity-weighted galaxy counts 
and uses only the masses of the haloes within a sphere. § [4] 
contains the main results of this paper, where we quantify 
how the merger rate is amplified in denser regions and sup- 
pressed in voids for redshifts 2: = to 2 over three decades 
of halo mass (10^^ — 5 x 10^'' ^^0). A simple power-law fit- 
ting function for this environmental dependence is proposed, 
which can be used in combination with the fit for the global 
rate presented in FM08. We also show that the progenitor 
(or conditional) mass function has a similar environmental 
trend as the merger rate. Even though this is expected given 
that the two quantities are closely related, this result demon- 
strates directly that the excursion set model is incomplete. 
In § 5, we present statistics of halo fragmentations, compare 
five algorithms for handling these events, and illustrate the 
robustness of the results reported in § 4. The Appendix pro- 
vides a discussion of the self-similarity of the merger rate and 
its environmental dependence in the context of the choice of 
mass and environment variables used in the fitting formula. 



HALOES IN THE MILLENNIUM 
SIMULATION 



The Millennium simulation (Springel e t al.|2005| ) follows the 
evolution of roughly 2 x 10^ dark matter haloes from redshift 
z = 127 to 2 = in a 500ft^^ Mpc box using 2160^ particles 
of mass 1.2xlO^M0 (aU masses quoted in this paper include 
the factor of h~^). It assumes a ACDM model with Q,m ~ 
0.25, fit = 0.045, JIa = 0.75, h — 0.73 and an initial power- 
law distribution of density perturbations with index n = 1 
and normalisation = 0.9. 



A friends-of-friends (FOF) group finder (Davis et al 



19851 with a linking length of 6 = 0.2 is used to identify 
haloes in the simulation. Each FOF halo (henceforth halo) 
thus identified is further broken into constituent subhaloes 
(each with at least 20 particles or 2.35 x 1O"'M0) by the 
SUBFIND algorithm which identifies gravitationally bound 
substructures within the host FOF halo (for more on SUB- 
FIND, see |Springel et al.|2001| |. 

The subhaloes are connected across the 64 available red- 
shift outputs to form a subhalo merger tree. Mergers are 
complicated processes and the particles in a given subhalo 
will not necessarily end up in a single subhalo in the sub- 
sequent output. As such, a subhalo is chosen to be the de- 
scendent of a progenitor subhalo at an earlier output if it 
hosts the largest number of bound particles in the progeni- 
tor subhalo. The resulting merger tree of the subhaloes can 
be used to construct the merger tree of the FOF haloes, al- 
though we have discussed at length in FM08 that this con- 
struction is non-trivial due to the fragmentation of FOF 
haloes. Our main results reported in Sec. 4 use the stitching 
tree of FM08. Since fragmentation occurs more frequently 
in denser environments, we provide a detailed comparison 
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Mass Percentile 



Redshift 


Quantity 


0-40% 


40-70% 


70-90% 


90-99% 


99-100% 


2 = 0(Az = 0.06) 


Number of haloes 
Mass bins (IO^^Mq) 
1/ bins 


192038 
1.2 - 2.1 
0.75-0.81 


144028 
2.1 - 4.5 
0.81-0.92 


96019 
4.5 - 14 
0.92-1.11 


43208 
14- 110 
1.11-1.63 


4800 
> 110 
1.63-4.30 


2 = 0.51 (Az = 0.06) 


Number of haloes 
Mass bins (IO^^Mq) 
1/ bins 


188258 
1.2 - 2.0 
0.95-1.03 


141194 
2.0 - 4.1 
1.03-1.15 


94129 
4.1 - 12 
1.15-1.37 


42358 
12 - 74 
1.37-1.93 


4706 
> 74 
1.93-4.66 


2 = 1.08 {A2 = 0.09) 


Number of haloes 
Mass bins (IO^^Mq) 
1/ bins 


172568 
1.2 - 1.9 
1.22-1.32 


129426 
1.9 - 3.7 
1.32-1.46 


86284 
3.7-9.5 
1.46-1.70 


38827 
9.5 - 48 
1.70-2.27 


4314 
> 48 
2.27-4.73 


2 = 2.07 (Az = 0.17) 


Number of haloes 
Mass bins (IO^'^Mq) 
V bins 


116830 
1.2 - 1.8 
1.74-1.86 


87622 
1.8 - 3.0 
1.86-2.01 


58415 
3.0 - 6.5 
2.01-2.27 


26286 
6.5 - 24 
2.27-2.85 


2920 
> 24 
2.85-5.20 



Table 1. Halo mass bins and number statistics at redshifts z = 0,0.51,1.08 and 2.07 from the Millennium simulation used in this 
paper. The bins are computed assuming fixed mass-percentile bins (header row). Listed are the number of haloes in each bin, and the 
corresponding mass and v boundaries for each percentile bin. The highest 1% mass bins extend out to 5.2 X 10^^,3 X 10^^,1.3 X 10^^, 
and 4.4 X IO^^Mq for 2 = 0, 0.51, 1.08 and 2.07 respectively. 



in Sec. 5 between stitching and four alternative algorithms 
to test the robustness of our results. 

The Millennium database provides a number of mass 
measurements for each identified FOF halo. We use the to- 
tal mass of the particles connected to an FOF by the group 



finder. Tinker et al. ( 2008 1 argue that spherical overdensity 
measures of mass are more closely linked to cluster observ- 
ables than FOF measures and, therefore, are to be preferred. 
We have found, however, that the FOF mass definition is 
more robust in the context of merging haloes than definitions 
that make assumptions about halo geometry and virializa- 
tion (for the simple reason that merging haloes are typically 
not virialized at the simulation outputs immediately preced- 
ing and following a merger event; see also White 200T] ). 

We study the dependence of halo growth on halo envi- 
ronment in a variety of halo mass bins at different redshifts. 
Since the most massive haloes at 2 = are more massive 
than the most massive haloes at higher redshifts, we use 
mass bins with boundaries that vary with redshift such that 
each mass bin contains a fixed percentage of haloes. Ta- 
ble [1] lists the five percentile bins used in our study, and 
the corresponding number of haloes and halo masses at 
z — 0.51, 1.08, and 2.07. We note that even in the highest 
1% mass bin, there are 4000 to 5000 cluster-mass haloes at 
z < 1 available for this study. 

Table[l]also lists the range of v for each mass bin, where 
f = Sc{z)/a{M) is often used as a mass variable for compar- 
ing haloes over different redshifts. Here a{AI) is the variance 
of the linear density perturbations and 5c{z) is the critical 
overdensity at redshift z, where 5c{z) — 1.686/D(2:) and 
D{z) is the linear growth function of density perturbations 
in ACDM. A comparison of the FOF mass versus v as the 
mass variable is provided in the Appendix. 

Our notation is as follows. When computing the merger 
rate, we refer to the haloes at the lower redshift as the de- 
scendants and label their masses by Mq. The progenitors 
of a given descendant halo at a (slightly) higher redshift 
are labelled Mi, M2, Mi..., where Mi > M2 ^ M3... by 
our convention. The mass ratio of the progenitors is de- 
fined as 5 = Mi^2/Mi. In this paper we find that there 



are sufficient halo statistics from the Millennium simulation 
for studying the environmental dependence of the descen- 
dant haloes over a range of redshifts, and shall present re- 
sults at zq — 0,0.51, 1.08, and 2.07 and their progenitors at 
zi = zo + Az, where Az = 0.06, 0.06, 0.09, and 0.17, respec- 
tively. 



3 MEASURING HALO ENVIRONMENT 

In this paper we quantify a halo's local environment using 
the local mass density centred at the halo. In this section 
we examine four definitions of density. Three of them are 
computed using the dark matter particles in a sphere of 
radius R centred at a halo, either with or without the central 
region carved out (see Sec 3.1-3.3). The fourth definition is 
computed using the masses of only the haloes rather than all 
the dark matter (Sec. 3.4). This last environmental measure 
based on mass-weighted halo counts has the advantage that 
it can be linked to observables such as luminosity-weighted 
galaxy counts. 



3.1 Definitions of Environment 

Only a few studies of halo environment have used local over- 



densities as measures of environment ( Lemson & Kauffmann 
[Harker et ah] [20061 |Wang et al.||2007[ |Hahn et al 



1999 



2008 1. By contrast, many studies have used the halo bias as a 
proxy for halo environment, which is obtained by taking the 
ratio of the halo-halo (or halo-mass) two-point correlation 
function to the underlying dark matter two-point correla- 



tion function (e.g.. 


Gottlober et al. 


2002 |Sheth & Tormen| 


2004||Gao et al. "2005; 'Marker et al. 


20061 IJing et al.||2007| 


Wechsler et al.,,2006, Gao & White 


2007 


1. Typically, these 



studies explore the dependence of bias on a variety of tracers 
of the halo growth history such as formation redshift, con- 
centration, and number of major mergers. This technique 
has yielded clear signs of environmental dependence, partic- 
ularly when combined with the marked correlation function 
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Figure 1. Scatter plots of halo mass vs three measures of halo environment for all FOF haloes above 1.2 X IO^^Mq (1000 particles 
or more) in the 2 = Millennium simulation output. The colour scale indicates the number of haloes present in each ((5, M) grid cell 
normalised by the bin size, and the contours are drawn at the 1, 100, and 10* bin levels (decreasing line width). The left panel uses 
1 + ^7, the density in a sphere of radius 7h~^ Mpc centred at each halo. The black line is 1 + ^7 = M /V-j / pm (see text). The middle 
panel shows 1 + 67-2, the density in a shell between 2h~^ and 7h~^ Mpc. The right panel uses 1 + (57_fof t>y subtracting the halo mass 
from St. At the high mass end, the halo itself is the main contribution to <57 and <57_2, leading to the tight correlation between & and 
M in the upper right region in the left and middle panels. The right panels shows that this correlation is largely removed when <57_pop 
is used, which subtracts out the FOF mass of the central halo. The variable S^—yoy is therefore a more independent measure of the 
immediate environment outside of the haloes. 



Statistical test ([Gottlober et al.|2002||Sheth fc Tormen|2004 
[Harker et al.|2006| [Wechsler et al.|2006| ). 

The connection between a halo's local density and the 
halo bias, however, is not entirely straightforward. The two 
quantities are certainly correlated, e.g., the two-point corre- 
lation function of objects in denser regions is typically higher 
than that in less dense regions ( [Abbas fc Sheth|2005| ). How- 
ever, the local density is a simple quantity that can be com- 
puted for each halo, whereas the bias is a statistical measure 
of clustering strength computed by averaging over a large 
number of pairs of haloes and particles over a range of pair 
separations. 

The rich statistics of the Millennium simulation over 
large dynamic ranges in both mass and redshift make it 
possible to use the more intuitive local density as a measure 
of environment. 

To compute the local overdensity in a halo's neighbour- 
hood, we centre either a sphere or shell on the halo at spatial 
coordinates x and define the halo's environment by 



5ij(x) 

for a sphere of radius R, or 



m-Rf 



(1) 



(2) 



for a shell of inner and outer radii Ri and Ro- Here pm is the 
mean matter density in the simulation box, and Pb(x) is the 
mean density of a sphere of radius R centred at x. We also 
propose an environmental measure, Sr-fqf, computed by 
subtracting out the FOF mass M of the central halo within 
a sphere of radius R: 



Sr- 



M 
Vrp„ 



(3) 



where Vr is the volume of a sphere of radius R. Note that 
unlike the shell measure, this measure makes no assumption 
about the central halo's shape. 



To compute p_b(x), one would need all the particle po- 
sitions from the Millennium simulation, which are not avail- 
able on the online public database. The database, however, 
does provide the density on a 256^ cubic grid (with a grid 
spacing of 1.95 Mpc) computed from the dark matter 
particles in the simulation using the Cloud-in-Cell (CIC) 
interpolation scheme. We use this data to sum up the con- 
tributions within the sphere centred at x to evaluate p_r(x). 
We note that the grid in the database is indexed using a 
Peano-Hilbert space filling curve, which we have mapped to 
spatial coordinates in order to compute pr. 

It is important to choose an appropriate radius R (or 
Ro — Ri) in equations ([lt-(|3l) when computing halo envi- 



ronment 

55-2 (in units of h 



Lemson fc Kauffmann ( 1999[ ) used both 5io and 



Mpc) but failed to detect any envi- 
ronmental dependence in the formation redsh ift for haloes 
with masses between 2 x 10^'^ and lO^*/i~^M0. Harker et al. 
( 2006) , following Lemson & Kauffmann (19991, used S5-2 



and did detect environmental dependence in Millennium for 
haloes with masses between 2 x 10^^ and lO^/i-^Mo. simi- 
larly, Hahn et al. ( 2008 1 used and 65-2 and found environ- 



mental dependence for haloes with mass es be tween 2 x 10 
and 1.6 x lO"/i"^Af0. We will show in ^S^sjthat R = 7/i"^ 
Mpc is an adequate choice that effectively characterises the 
environments of massive halos. 



3.2 Disentangling Environment and Mass 

Since the goal of this paper is to quantify the dependence of 
merger rates on halo environment as well as mass, it is es- 
sential for us to first examine the extent to which these two 
variables are independent measures of halo properties. This 
is particularly relevant considering that our measure of en- 
vironment is based on the local mass density. We note that 
in the literature mass is often used loosely to refer to envi- 
ronment, e.g., clusters are considered denser environments 
than galaxies. This interpretation is valid for galaxy counts. 
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We are concerned with FOF haloes (and not subhaloes or 
galaxies) here, however. As we will show, haloes of all masses 
can reside in a wide range of overdensities. 

To study the relation between halo mass and environ- 
ment, we present a scatter plot of the mass of every FOF 
halo (above 1000 particles M > 1.2 x lO^^M©) at z = in 
the Millennium simulation versus its local 1 + 5 in Fig. [l] 
Three definitions of local density are shown for comparison: 
all mass within a 7h~^ Mpc sphere (^7; left panel), all mass 
within 7h~^ Mpc excluding the central 2h~^ Mpc ((57_2; 
middle panel), and all mass within 7h~^ Mpc excluding the 
central FOF mass (57_fof; right panel). 

Fig.[l]shows that galaxy-size haloes (~ 10^^ Mq) reside 
in a wide range of environmental densities from extreme 
underdense regions of 5 ~ —0.8 to regions with 5 > 20. 
The three panels show similar distributions of S for these 
low mass haloes regardless of the definition of 5 used. The 
high mass haloes, on the other hand, have very different 
distributions of 5. The rich statistics of the Millennium sim- 
ulation allow us to study halo mass out to 2 x W^^Mq, 
an order of magnitude higher than in previous studies. At 
5 X 1O^*M0 and above, the spherical and shell measures of 
overdensity, 87 and ^7-2, are seen to be tightly correlated 
with the halo mass (left and middle panels). In addition, 
all the points lie close to the line that represents the den- 
sity in a 7h~^ Mpc sphere computed from the FOF halo 
mass alone, that is, M /V? / pm, where V? is the volume of 
a sphere of radius 7h~^ Mpc. This trend clearly indicates 
that the central haloes are dominating the local overdensity 
at masses above ~ 5 x 10^*Mq, and both and ^7-2 are 
tracing the central halo mass rather than the overdensities 
in the neighbourhood outside the virial radius of the halo. 
Even though 67-2 subtracts out the central 2h~^ Mpc re- 
gions, the tight residual correlation seen in the middle panel 
suggests that this quantity does not cleanly remove the con- 
tribution made by the central halo, probably because these 
massive haloes extend well beyond 2h~^ Mpc. We have also 
tested S5-2, the measure of environmen t used in|Lemson fc| 
Kauffmann] ( |1999[ ); [Harker et aL| ( |2006[ ); |Hahn et al.| ( |2007f 7 
and found nearly identical results as ^7-2. This correlation 
may not be problematic for the results reported in these ear- 
lier papers, however, as these studies did not report results 
beyond ~ IO^Mq. 

The right panel of Fig. [l] shows that our third envi- 
ronmental variable, 57_fof, in equation ([sjl is capable of 
disentangling the tight correlation between halo mass and 
density seen for S7 and 87-2; S7-P0F is therefore a more ro- 
bust measure of the environment outside of a halo's virial 
radius. It should, however, be kept in mind that haloes of 
different masses residing in the same 7 Mpc region will have 
the same 87 but different (57_fof. A cluster-sized halo, for 
instance, will have a smaller value of (57_fof than a neigh- 
bouring galaxy-sized halo, and the difference between the 
two values of S7-F0F will be the difference between the mass 
of the cluster and the galaxy (appropriately normalised). 
This caveat should be considered when interpreting values 
of (57_FOF across different mass bins. The spherical measure 
67, on the other hand, is simpler in this context. For this 
reason, we will report results using both St-fqf and S7 be- 
low. 



3.3 S Distributions 

To gain further insight into the properties of the environ- 
mental measure Sr-fof of equation (|3|, we plot in the left 
panels of Fig. [2] the distribution of 1 4- 5r-fof centred at 
each halo for all haloes in the Millennium database with 
more than 1000 particles (M > 1.2 x lO^^M©) at four red- 
shifts z = 0,0.51,1.08, and 2.07 (top to bottom). Within 
each panel, four choices of radii, R — 3,5, 7, and 9h~^ Mpc, 
are shown for comparison (thin black, thin dark grey, thick 
black, and thin light grey). A comparison of the four left 
panels shows that the width of the 1 -I- (5h-fof distribution 
becomes broader towards lower redshifts. This is a natu- 
ral consequence of gravitational instability: denser regions 
become denser and vice versa as the universe evolves. We 
will explore the implications of this effect further in the ap- 
pendix. 

At a given redshift, as expected for a ACDM model, the 
overdensities computed using a larger smoothing radius R 
are generally smaller than those computed using a smaller 
R. In the voids, the distribution of 1 -I- Js-fof is seen to 
have a low 5 tail that extends down to unphysical (nega- 
tive) 1 + 5; a, faint remnant of this tail is also visible in 
1 -I- (Js-fof at z = 0. This tail is due to a number of cluster- 
size haloes whose FOF member particles extend beyond 3 
to 5 Mpc. We note that the virial radii of even the most 
massive halos (lO^^A/0) do not extend beyond 3 Mpc; how- 
ever we have found that the distance between the centre of 
an FOF's most massive subhalo and the furthest subhalo 
associated with said FOF can extend beyond 5 Mpc even 
for halos with masses of a few xIO^Mq. We therefore use 
1 -I- 57_FOF throughout this paper in order to better sample 
the environment surrounding these large haloes. 

We break down the haloes represented by the thick 
black 1 + 57-FOF curve in the left panels of Fig. [2] into dif- 
ferent mass bins and plot their (5-distributions using dot- 
ted colour histograms. Less massive haloes cover a broader 
range of 1 -I- (57_pof than more massive haloes, and their 
distribution peaks at a lower value of (57_pof. We note that 
even though the histograms for both the lower mass haloes 
and the total distribution at z — are peaked at a slightly 
negative value of S, the mean value is in fact positive, e.g., 
< 57 >= 0.864 and < (57_pof >= 0.849 for the lowest mass 
bin, and < 87 >= 1.04 and < 57~fof >= 0.956 for all 
the halos. The value of < 57_fof > is only slightly smaller 
than < 87 > because subtracting the FOF mass of the low 
mass haloes (which dominate the total distribution) makes 
little difference when 5 is averaged over a sphere of radius as 
large as 7h~^ Mpc. The mean of S is not zero here because 
the overdensities are not randomly sampled but are instead 
centred on haloes. 

The right panels in Fig.[2]are scatter plots of each halo's 
local density 1 + (57_pof versus its mass a.t z — 0, 0.51, 1.08, 
and 2.07 (top to bottom). (The top panel is a repeat of the 
right panel of Fig. [l]) The mass bins based on percentiles 
from Table [1] are marked by the horizontal lines. At high z 
the haloes cover a narrower range in both S and M, but the 
tight correlation seen in Fig.[l]between the mass of the mas- 
sive haloes and their local densities ^7 and 87-2 is removed 
at all redshifts when J7_fof is used. 
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Figure 2. Left panels: Distribution of tlie environmental variable 1 + Sr-fof defined in equation (jSjl at four redshifts z = 0,0.51, 1.08 
and 2.07 (top to bottom) computed from all haloes with M > 1.2 X 10^^ (i.e. above 1000 particles) in the Millennium simulation. The 
broadening of the distribution with decreasing z illustrates the effect of gravitational instability. Within each panel, the four grey-scale 
histograms compare four smoothing radii R in Mpc: 3 (thin black), 5 (thin dark grey), 7 (thick black) and 9 (thin light grey); the 
five coloured dotted histograms compare the separate contributions to the S^-poF distribution from haloes of different mass percentile 
bins: top 1% (red), 90-99% (olive), 70-90% (green) 40-70% (cyan), and bottom 40% (blue). Right panels: Similar scatter plot as the right 
panel of Fig. [T] but at four redshifts. The horizontal dotted lines mark the five mass percentile bins used in the left panels and listed in 
Table [U 
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Figure 3. Scatter plot of two density variables; 1 + (see 
Fig. 1) computed from all dark matter particles centred within 
a 7h~^ Mpc sphere of each halo, and the mass- weighted halo 
counts 1 + computed by including only masses in haloes in 
the same sphere (down to halo mass of 1.2 X IO^^Mq). The con- 
tours are plotted at probability values of 10~^ (black), 0.01, 0.1 
and 1 (white) ; the thick contours are for simulation data, the thin 
contours are from the fit in eq. The grey dotted line is for 
57 = 5f . 



3.4 Computing S via Halo Counts 

The local environmental measure S7 is convenient from a 
theoretical standpoint but is not easy to measure observa- 
tionally as it demands accurate knowledge of the background 
dark matter distribution within a large (7 Mpc) radius of the 
halo in question. Here we consider a more observer- friendly 
quantity based on the mass-weighted halo counts (above a 
certain mass threshold): 



l + <5f (x) 



(4) 



where the sum is over all halos within a 7h~^ Mpc sphere 
centred at x above some minimum mass (we use 40 particles, 
or 4.7 X 10^° M0), and W is the volume of a sphere of radius 
7h-^ Mpc. For a given halo in the Millennium simulation, 
we compute this quantity by summing over all haloes whose 
centres lie within the 7h~^ Mpc sphere centred on the halo 
in question. We do not account for the fact that halos near 
the boundary may only strictly contribute a fraction of their 
mass to the 7h~^ Mpc sphere. 

The resulting mass- weighed halo counts 1 + is plot- 
ted against 1 + 67 computed from the CIC density grid in 
Fig. |3] The 2d-histo gram is normalised to have unit area 
and can be thought of as a bivariate probably distribution. 
We note that, while Sr is generally greater than 3^' as ex- 
pected, there are regions with 5^ > 3?, particularly in dense 
environments. This is due to the fact that a halo's entire 



mass contributes to 3^ if its centre lies within the 7h~^ 
Mpc sphere in question. 

We approximate the distribution with a two- 
dimensional log-normal distribution. Since the variables are 
correlated, the fitting form has five parameters and is given 
by 



dP 



exp 



ln(xi 



ln(a;2) 



2Gi 



2ol 



(5) 



where In(xi) and ln(a;2) are uncorrelated variables that are 
simply linear combinations of ln(l-f57) and \a{\^37^) given 
by 



ln(a;i) 




cos Q 


sinS 




ln(l -h (57) - Ail 


ln(x2) 




— sin Q 


cos d 




ln(l + 5f ) - /i2 



(6) 

Here /xi and /i2 denote the mean values of ln(l -I- ^7) and 
3^^) respectively, Q is an angle that quantifies the cor- 
relation between the two (5s, and a\ and (72 are the standard 
deviations along the major and minor axes defined by Q. 
The best fit values for these five parameters are /ii = 0.210, 
[1-2. = -0.549, (Ti = 1.03, (72 = 0.141, e = 0.943. The thin 
contours in Fig. [S] represent the resulting fit. 

We have also computed a simpler power-law fit that can 
be used to approximate the mean relation between the two 
densities: 



ln(l -f 5f ) = 1.28 ln(l + 57) - 0.865 . 



(7) 



Both fitting forms can be used to convert back and forth 
between 3^ and 3^^ . 



4 ENVIRONMENTAL DEPENDENCE 
4.1 Halo Merger Rate 

In FM08 we defined and computed the merger rate Bjn 
as a function of progenitor mass ratio f = Mi/M\ (with 
i ^ 2), descendant mass Mo, and redshift z. The rate B/n 
is dimensionless and measures the mean number of mergers 
per halo per redshift interval per mass ratio. We found that 
in these units, the merger rate has a remarkably simple form 
and depends only weakly on mass and redshift. We proposed 
the fitting form 

f)"^=-[(l)"l(§ 



B{Mo,i,z) 



= A 



n{Mo,z 
where (Q,/3,7,r;) = (0.083 



(8) 



-2.01,0.409,0.371), A = 0.0289, 
i = 0.098, M = 1.2 X IO^^Mq, and 5c{z) oc 1/D{z) is 
the standard density threshold normalised to 5c ~ 1.686 
at 2: = 0, with D{z) being the linear growth factor. This fit 
is accurate to 10-20% over the mass range lO'^^ — IO^^Mq 
and redshift range z < 6. 

The merger rates B/n at z — for the five descendant 
halo mass bins in Table [l] are reproduced for reference in the 
top left panel of Fig. [4] As shown in FM08 and indicated by 
equation the merger rate is approximately a power law 
in ^ in the minor merger regime and has a slight upturn in 
the major merger regime > 0.2). All five curves are nearly 
on top of one another, reflecting the very mild mass depen- 
dence (q ~ 0.1) in equation ([8|. We compute each curve 
by first selecting the descendant haloes in a given mass bin 
and computing the mass ratios ^ for the progenitors of these 
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Figure 4. Halo merger rate and its environmental dependence on the local overdensity (57_poF measured in a 7h~^ Mpc sphere excluding 
the central FOF halo mass. Top left panel: The global mean merger rate B/n (in units of mergers per descendant halo per unit redshift 
per J bin) as a function of the progenitor mass ratio ^ for descendant haloes in five mass percentile bins (see Table Q . The results are 
computed using the 2 = and 0.06 outputs from the Millennium simulation. The higher mass curves extend down to lower ^ because we 
have chosen a fixed minimal progenitor mass (40 particles) for all descendants. Other five panels: The ratio of the merger rate of haloes 
in a given environmental bin B/n[S\ to the global mean B/n as a function of ^. Each panel is for a mass bin shown in the upper left 
panel. Within each panel, different colours show different 1 + 57_poF bins (red for the densest regions; blue for the void regions), and 
the bands indicate the size of the Poisson errors. Note that the lower panels for the higher mass haloes have fewer & curves since 5-j—yO¥ 
for these haloes spans a narrower range. This figure clearly shows that the merger rate is higher in dense regions and lower in voids for 
all halo masses, and the boost or reduction factor is nearly independent of the progenitor mass ratios 



haloes. We then compute B/n by counting B, the number 
of progenitors [Mi with i > 2) that He in a given mass ratio 
bin (^), and dividing by n, the total number of descendants 
in the mass bin in question. See FM08 for further details of 
this procedure and discussions of the results. 

Equation ([8| gives the global mean merger rate aver- 
aged over all halo environments. To investigate the corre- 
lation of B/n with environment, we divide each mass bin 
shown in the top left panel of Fig. [4] into five environmental 
bins and compute B /n[S\ using B and n in the given (Jt-fof 
bin. The remaining five panels in Fig. |4] show our results for 
the ratios of B/n[S\ to the global mean B/n, as a function 
of f , for each of the five descendant mass bins. Within each 
panel, the different curves are for different St-fof bins for 
which there are sufficient halo statistics. 

For descendant haloes of mass 10^^ to lO^M©, Fig. [i] 
shows a strong environmental effect with a positive correla- 
tion between merger rates and local density: haloes in the 
densest regions (1 -I- 5i-fof > 7; red curves) experience 1.5 
to 2 times more mergers than the average, while haloes in 
underdense regions (1 -I- ^t-fof < 0.7; blue curves) experi- 
ence fewer mergers (by a factor of 0.7 to 0.8) than average. 
For group and cluster scale haloes (10^"^ to 5 x lO^^M©) in 
the bottom panels, only three curves are shown for the mid- 



dle three (Jt-fof bins because massive haloes span a smaller 
range of (Jv-fof, as shown in Figs. 1 and 2. For each 5t-fof 
bin, the value of the merger rate ratio is quite similar for 
all five panels, indicating that the environmental effect, as 
measured by Sj-fqf, depends very weakly on halo mass. We 
will quantify this statement using a fitting formula below. 

Fig. [5] presents the same information as Fig. |4] but at 
higher redshifts {z — 0.51, 1.08, and 2.07 from top to bot- 
tom panels). The 5-bins for which the curves are too noisy 
are excluded. Since the distribution of 1 -f S^-fof narrows 
with increasing 2, the S7-F0F bins span a smaller range at 
z = 2 than at 2 = 0. Nonetheless, we see that the environ- 
mental dependence observed at z = persists out to 2: = 2. 
Moreover, haloes with similar 1 -|- 57_fof experience similar 
amplifications or reductions in the merger rate regardless of 
mass and redshift. 

An additional feature to note in Figs. |4] and [5] is that 
the curves are horizontal: environmental effect is therefore 
largely independent of the mass ratio ^; that is, major and 
minor merger rates are boosted or dampened by a halo's 
environment by a similar factor. We can therefore integrate 
over the mass ratio parameter without diluting the environ- 
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Figure 5. Same as the merger rate ratio plots in Fig. [4]except at higher redshifts: z = 0.51, 1.08,2.07 (from top to bottom). The five 
columns correspond to the five mass percentile bins (see Table [l]). Within each panel, the coloured curves show different 1 + 57_pop 
bins (red for the densest regions; blue for the void regions), and the bands indicate the size of the Poisson errors. Note that since the 
distribution of (57_fof evolves with z, the corresponding 5 bins (labelled in the leftmost columns) for the five coloured curves change 
with redshift. This figure shows that the positive correlation of the merger rate with <57_poF persists out to ^ ~ 2. 



mental effect: 



dz 



'-(M,z) 



n{M, z) 



(9) 



where dNyncr^c/ dz is the mean merger rate per unit redshift 
per descendant halo with progenitor mass ratio above ^min- 



The value of dNmcrgc/ dz clearly depends on ^min and 
is larger when more minor mergers are included (see, e.g., 
Figs. 7 and 8 of FM08). For a fixed resolution mass (our 
choice is 40 particles or more for progenitor haloes), ^min 
extends down to lower values for higher mass descendants. 
For a fair comparison across halo mass bins, one should in 
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Figure 6. Dependence of the mean merger rate dNjncigc / dz (= J B/ndi;) on environmental variables 1 + 57_poF (top figure) and l + Sr 
(bottom figure) at four redshifts z = 0, 0.51, 1.08, and 2.07 (left to right). Within each figure, the top panel shows the ratio of the mean 
merger rate dNmcrgc / dz[5] for haloes in a given environment to the global mean merger rate dNmcrgc/dz. The bottom panel plots the 
ratio of the simulation results to the fits, showing that eq. is generally accurate to within 10% (indicated by the dotted horizontal 
line). The colours correspond to the five mass percentile bins in Table[T](red for the highest and blue for the lowest mass bin); the bands 
correspond to Poisson errors. This figure shows that the positive correlation of the merger rate with environmental density is present at 
all mass and redshift ranges probed. 



principle use a fixed ^min for all mass bins at the expense of 
throwing out resolved progenitors for high mass descendant 
haloes. Since we plot ratios of the merger rates, however, 
the fact that more massive haloes are better resolved and 
have higher dNmeige/dz is normalised out. It is therefore 
possible to make a fair comparison across mass bins without 
throwing out any resolved progenitors. 

Fig.[6]shows the ratio of the merger rate dNmcrgc/dz as 
a function of 1 + 5 at four redshifts {z = 0,0.51,1.08,2.07 
from left to right). For comparison, the results for two envi- 
ronmental measures are included: ^t-fof (top figure) and 5r 
(bottom figure). Within each figure, the upper panel shows 
the simulation data and the lower panel compares the data 



to the analytic fitting formula discussed below. The five 
curves in each panel are for the five mass bins listed in Ta- 
ble [l] This fi gure shows the same trend as Fig. [4] haloes in 
the densest regions at 2 = experience up to ~ 1.5 times as 
many mergers as the average halo, whereas the merger rate 
in the voids is 20 to 30% below the global average. 

To quantify the dependence of the merger rate on 5, we 
introduce 

^^^{S,M,z)^^^^{M,z)xf{5,M,z), (10) 

where we have made use of the fact that the environmental 
dependence is independent of ^ to define /, and dN / dz{M , z) 
is the global merger rate from FM08. We provide two fit- 
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Figure 7. Same as Fig. |4] except the quantity shown is the progenitor (or conditional) mass function (l>(AIi, zi\Mo, zq) instead of the 
merger rate B/n. The results are computed using the (^Oi^i) = (0,0.06) outputs from the Millennium simulation. Similar to Fig. |4] we 
find that descendant haloes in the densest regions (red curves) have significantly more progenitors (by a factor of 1.5 to 2) than the 
global distribution of progenitors, whereas those in the voids (blue curves) have ~ 20% to 30% fewer progenitors. 



ting forms for / using 5 = i^t-fof and 5-^, respectively. We 
find that a simple power-law and redshift-independent form 
works well: 



of 57_FOF, as the environment variable, the only change in 



/(<57-FOF, M) = 0.963 (1 + 57-fof) 



M 



f{5-j,M) = 0.968(1 + ^7) 



M 



1012 Afg 



1012 M© 

-0.0252 



(11) 



The fits are performed using data from all four redshifts si- 
multaneously (over much finer mass bins than those shown 
in Fig. [6]). The reduced for the two fits is 0.95 and 1.08, 
respectively. Errors are computed assuming Poisson statis- 
tics and are represented by the filled regions in Fig. [6] The 
resulting fit is shown as dashed curves in the upper row of 
each figure, and the ratio of the simulation data to the fits 
is shown in the lower rows. The fits are seen to be accurate 
to within 10% for a wide range of 5 and M, except for low 
mass haloes with (5 > 5 at z = and 0.51, where the rates 
steepen suddenly. 

As we will discuss in § |5] our extensive tests using vari- 
ous algorithms suggest that the merger rate in this particular 
parameter range (i.e. low mass, low z, high density) depends 
sensitively on the post-processing algorithm used to handle 
fragmentations in the merger tree, and variations of order 
20% or more among different algorithms are observed. We 
therefore do not attempt to use a fitting form more compli- 



cated than equation (111 to get a better fit in this uncertain 
regime. 

It is interesting to note that when S-j is used, instead 



the fit in equation (111 is a stronger dependence on halo 
mass. This trend makes sense since the difference between 
^7 and (57-FOF is 5-j - &7-fof = M/{V7pm) (see eq. [s]). 
This difference is negligible for galaxy-scale haloes (e.g. 5-j — 
57_FOF ~ 0.01 for IO^^Mq) but becomes larger for more 
massive haloes, reaching ^7 — Sj-yof ~ 10 at M ~ 10^'^Mq. 
The five curves for the five mass bins at a given z in Fig.[6]are 
therefore more spread out when ^7 is used as the variable, 
resulting in a stronger mass dependence. 

We have chosen to use halo mass and local density as 



variables in equation (111. It is interesting to ask if other 



choices of variables may lead to a more accurate fit across 
the wide ranges of halo masses, densities, and redshifts 
shown in Fig. [6] For instance, the variance of the linear 
density perturbation a{M) and the scaled density threshold 
i^(M, z) — 5c/cr{M)D{z) are commonly used to characterise 
mass and redshift dependence of halo properties (e.g., the 
mass function). We test these variables and describe the re- 
sults in appendix]^ Our conclusion is that these alternative 
variables do not perform any better, and the [M, 1-|-57_fof] 
pair shows the least systematic variation with redshift. 

In summary, the simple parametrisation of the envi- 
ronmental dependence of the merger rate given by equa- 



tions (10 1 and (111 can be used along with the fit for the 
global merger rate B/n in equation ([sjl to compute the 
merger rate in different environments at a variety of red- 
shifts. 
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4.2 Progenitor Mass Function 

For completeness and ease of comparison with analytic mod- 
els, we present liere the results for the environmental de- 
pendence of the conditional (or progenitor) mass function 
(/)(Mi, zi\Mo, zo). This function gives the mean distribution 
of the progenitor masses Af; at redshift Zi for a descendant 
halo of mass Mq at redshift zq. It is the key ingredient for the 
construction of Monte Carlo merger trees in the Extended 
Press-Schechter model. 

The relation between (f) and the merger rate B/n is dis- 
cussed in Sec 3.3 of FM08. These two quantities are closely 
related but differ in two ways. First, is typically plotted 
vs Ml/Mo, while B/n is expressed in the mass ratio of the 
progenitors Mi/Mi{i ^ 2) and the descendant mass Mq. 
Second, the conditional mass function (j){Mi, zi\Mo, zq) in- 
cludes all progenitor halos at zi regardless of if a merger 
has occurred between zi and zo, whereas the merger rate 
includes only descendant haloes with more than one pro- 
genitor. When the lookback time zi — zo is small, a large 
fraction of haloes in fact have only one resolved progeni- 
tor typically with a mass Mi comparable to the descendant 
mass Mq. See the sharp rise in (f) near Mi/Mq = 1 in Fig.[7|. 
No such peak is present in the merger rate in Fig. |4] 

Fig. [T] shows that the progenitor mass function has a 
similar dependence on 1 -I- iSt-pop as the merger rate in 
Figs . |4)6"1 We have chosen to plot Fig. [7] in the same way as 
Fig. H] where the upper left panel shows the global progeni- 
tor mass function (f) aX z\ = 0.06 for five bins of descendant 
mass Mo at zq = 0, and the other five panels show how 
cf> for haloes in different (St-pop bins compare to the global 
mean (f>. We see that, like B/n, the progenitor mass function 
has a noted dependence on environment. For galaxy-size de- 
scendant haloes, those in the overdense regions have ~ 1.5 
times as many progenitor haloes as the mean, while those in 
the underdense regions have ~ 0.7 times as many progenitor 
haloes as the mean. 



5 ALTERNATIVE ALGORITHMS FOR 
POST-PROCESSING HALO 
FRAGMENTATIONS 

As we discussed in Sec. 2 and FM08, even though each sub- 
halo in the Millennium tree is, by construction, identified 
with a single descendant subhalo, the resulting FOF tree 
can contain fragmentation events in which an FOF halo is 
split into two (or more) descendant FOF haloes. This frag- 
mentation issue is not unique to the use of subhaloes in 
the Millennium simulation but occurs in merger trees in all 
prior studies that are typically constructed based on the 
FOF haloes rather than subhaloes. This problem arises be- 
cause particles in a progenitor halo (or subhalo) rarely end 
up in exactly one descendant halo; a decision must therefore 
be made to select a unique descendant and there is no unique 
way to do this. The standard procedure to assign progenitor 
and descendant FOF haloes is the same as that applied to 
the subhaloes in Millennium: the descendant halo is the halo 
that inherits the most number of bound particles of the pro- 
genitor. We call this algorithm snipping since it effectively 
cuts off the ancestral link between a progenitor halo and its 
subdominant descendant fragments, while leaving the halo 
masses unchanged (see Fig. [8|. 
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Split 






Figure 8. Example of a typical fragmentation event in the Millen- 
nium simulation. Black circles represent FOF haloes; white circles 
represent subhaloes. Circle radii scale with the logarithm of the 
(sub)halo mass. The left panel shows a fragmentation event oc- 
curring between z = 0.06 and 0.04 (red subhalo). The snip panel 
shows how the ancestral link between the fragmented halo and 
its progenitor is severed, producing a (blue) orphan halo. The 
stitch panel shows how the fragment is stitched back into the 
main branch at 2 = 0.04 (yellow subhalo). The split panel shows 
how the fragment's progenitor at 2 = 0.06 is split off from the 
FOF halo (green subhaloes). 



In FM08, we explored a new method stitching for han- 
dling these fragmentation events. In this method (which 
we call stitch-3 here), the fragmented haloes that remerge 
within 3 outputs after fragmentation occurs are stitched into 
a single FOF descendant; those that do not remerge within 
3 outputs are snipped and become orphan haloes. We com- 
pared the two methods and showed that snipping infiates 
the merger rates by up to 10% in the major merger regime 
and 25% in the minor merger regime (Fig. 9 of FM08). This 
is not surprising since bound subhaloes are often on eccen- 
tric orbits that extend out to 2 to 3 virial radii of the main 
halo (see, e.g., [Ludlow et al.pOOS] . The FOF finder can re- 
peatedly disassociate and associate these subhaloes, leading 
to spurious fragmentation and remerger events. 

In addition to the snipping and stitching algorithms, 
we examine a third method here that is complementary to 



stitching. We call this method splitting (see also Genel et al. 
(2008j )). Our motivation for introducing this algorithm is 
the fact that fragmentations can be the result of either false 
fragmentation at the lower zo, where physically bound sub- 
haloes are broken up, or false grouping at the earlier zi, 
where physically unbound subhaloes are falsely associated 
by the FOF finder. Even though multi-body subhalo encoun- 
ters may unbind a subhalo, our visual inspections of a num- 
ber of halo merger tracks indicate that such events are rare. 
Instead, most of the apparent fragmentations are due to the 
halo finder, which at an earlier output {zi) may group sub- 
haloes together, only to separate them at the next timestep 
{zq < zi). A decision needs to be made about whether the 
falsely separated haloes at zq should be put back together 
(i.e. stitching), or the falsely grouped halo at zi should be 
broken up (i.e. splitting). (Note: Snipping effectively does 
nothing.) 

Fig. [8] illustrates how each of the three algorithms - 
snip, stitch, and split - handles halo fragmentation. For com- 
pleteness, we also explore a variation of stitch (and split), in 
which the number of outputs used to make the decision is 
altered. Instead of stitch-3 (or split-3), which only stitches 
(or splits) fragmented haloes that remerge within 3 time 
outputs, we consider stitch-oo (or split-oo), which stitches 
(or splits) any fragmented haloes regardless of their future 
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Figure 9. Fragmentation statistics (fraction of haloes that fragmented) as a function of environmental density 1 + <57_fof s-t redshifts 
2 = 0.06,0.51,1.08,2.07 (left to right). Within each panel, three fragmentation mass ratios § are shown: major fragmentations with 
^ > 0.1 (dotted); those with ^ > 0.04 (dashed), and all fragmentations down to the resolution limit (40 particles; solid). The two colours 
are for different mass bins: to 40% (blue) and 90-99% mass bin (red); see Table [l] The red solid curve is significantly higher than the 
blue solid curve because the fragments of higher mass haloes are better resolved (i.e. ^rnin is smaller). Minor fragmentations are seen to 
dominate, while only ~ 1% of haloes suffer major fragmentations (^ > 0.1) in typical environments. 



( ( ( ( 
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e e ? ? 

Figure 10. Comparison of the global merger rate B/n (including all environment) vs progenitor mass ratio § computed from five 
fragmentation algorithms at z = 0.06,0.51, 1.08, and 2.07 (left to right). Since stitch-3 is the method used in FM08, we plot the ratio of 
B/n from the other four methods to B/n from stitch-3. Colours correspond to the mass bins in Table [l] (blue: lowest mass, red: highest 
mass). The shaded regions denote Poisson errors. 



(or past) history. An important distinction between stitch-3 
and stitch-cxj (and similarly for split-3 vs split-oo) is that the 
modifications to the haloes are confined to the three adja- 
cent outputs in stitch-3 and split-3; haloes along the merger 
tree outside of this time range are unaltered. The modifica- 
tions made in stitch-cxa and split-cxa however, propagate in- 
definitely either forward or backward along any tree branch 
where a fragmentation occurs. No algorithm is perfect, but 
any error made in stitch-oo and split-cxj will affect the en- 
tire branch of the tree that contains a fragmentation event. 
By contrast, errors made in stitch-3 and split-3 are confined 
to the redshift at which the fragmentation occurs. Stitch- 
oo and split-oo are therefore extreme algorithms, which are 
included here for comparison purposes only. 

Before comparing the algorithms, we first show the fre- 
quency of fragmentations in the Millennium FOF tree as a 
function of environment in Fig.|9] The fraction of haloes that 



experience fragmentations is seen to increase with (Jt-fof, 
differing by a factor of ~ 3 at low z and by a factor of ~ 5 
to 10 at z ~ 2. The fragments, however, are dominated by 
low-mass haloes: only ~ 1% of the haloes in typical densities 
have fragments of mass ratio ^ above 0.1, and this fraction 
is no larger than ~ 10% even in the densest regions. Most 
of the fragmentations are therefore minor. 



Fig. 10 compares the global mean merger rates B/n 



(i.e. including all environment) as a function of progenitor 
mass ratio ^ for the five algorithms (top to bottom) at four 
redshifts {z = 0.06,0.51, 1.08,2.07 from left to right). Since 
stitch-3 is the algorithm used in FM08, we show the ratio of 
each of the four alternative algorithms to stitch-3. Within 
each panel, the coloured curves show a variety of descendant 
mass bins (the bands show Poisson errors). As we have al- 
ready seen in FM08, snipping (first row in Fig. \10\ yields a 
higher merger rate (by ~ 20% at ^ < 0.01) due to the or- 
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1 — iSt-FOF 1 + (^T-FOF i + (St-FOF 1 + ^7-FOF 

Figure 11. Comparison of the environmental dependence of the merger rate computed from five algorithms at z = 0.06, 0.51, 1.08, 2.07 
(left to right). Similar to Fig.|6] we plot the ratio of the mean merger rate dTVmerge/c'zl'S] for haloes in a given 5j—fof bin to the global 
mean merger rate dNmeige / dz . The mass bins are shown in different colours (blue for the lowest and red for the highest bin in Tablejljl. 
Shaded regions indicate Poisson errors. This figure shows that stitch-3, snip, and split-3 have similar 5 dependence. Split-oo, however, 
reverses the <5-trend for low mass haloes, which we believe is an artefact of the propagation of fragmentations up the tree (see text). 



phaned haloes, resulting in a steeper power-law dependence 
^f^ (/? ~ -2.2) than stitch-3 (/? ~ -2). Stitch-oo (second 
row) , on the other hand, zips together all the fragments and 
reduces the number of minor mergers by as much as ~ 40% 
[P ~ —1.8) in comparison to stitch-3. Split-3 (third row) 
tends to raise the minor merger rate by up to ~ 15 — 20%. 
Split-oo (fourth row) has a feature in the low-^ merger rate 
that breaks the power-law behaviour seen in the other trees. 
This feature is redshift dependent and drives the merger rate 
lower than in stitch-3. 

We now examine how the environmental dependence 
of the merger rates is affected by the algorithm used for 
handling fragmentations. To do this, we integrate B/n over 
^ and show the total rate, dNn-icr^a / dz , as a function of 
1 + (Jt-fof for the five algorithms (top to bottom) at four 
redshifts in Fig. |11[ Similar to Fig. |6j the vertical axis 
shows the ratio of the merger rate in a 57_fof bin to the 
global rate, dNn-icr^c / dz[5] : dNmei^e/dz, computed with 
each method. This figure shows that the stitching (both 
stitch-3 and stitch-oo) and snipping algorithms produce 
very similar environmental dependence; though the extreme 
stitch-oo yields a mildly weaker S dependence. In contrast, 
split-oo shows a sudden reversal in the 5-dependence at 
z — 0.51, 1.08, 2.07 in the three lower mass bins, with haloes 
in the densest regions experiencing fewer mergers. Split-3, 
on the other hand, shows a positive (albeit weak) correlation 
of merger rate with S for all mass bins but the very lowest. 
We beUeve the difference between the split and stitch trees 
is due to an "unzipping" effect that is most pronounced in 
split-oo, in which splitting a fragmentation event at low z 
affects the entire branch above this redshift, resulting in the 



discrepantly low merger rates in high density regions seen 
in the last row of Fig. [TT] 

In summary, halo fragmentation is a generic feature of 
all merger trees. It occurs more frequently in dense regions 
than in voids, thereby prompting the detailed investigation 
in this section. Our tests of five algorithms show that the ma- 
jority of tree-processing methods (stitch-3, stitch-oo, snip, 
and, to some extent, split-3) give very similar environmen- 
tal dependence for the mean merger rate. In addition, the 
global merger rate (including all environment) is robust, dif- 
fering by less than 10% for major mergers and less than 20% 
even in the very minor merger regime (^ < 0.01) that is more 
prone to systematic effects. The split-oo algorithm, on the 
other hand, appears to suffer from "non-local" effects that 
have propagated up the merger tree from the fragmentation 
point. In particular, the merger rate is greatly reduced in 
high density regions when split-oo is used. 



6 CONCLUSIONS AND IMPLICATIONS 

We have used the dark matter haloes and merger trees con- 
structed from the Millennium simulation to quantify the de- 
pendence of halo merger rates on halo environment from 
redshift z = to 2. A number of local mass density parame- 
ters centred at the haloes, both including and excluding the 
central halo mass itself, are tested as measures of environ- 
ment. We have found that ^t-fof defined in equation ([3| is 
a robust measure of the surrounding environment outside of 
a halo's virial radius. It cleanly subtracts out the contribu- 
tions to the local density from the central halo and thereby 
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breaks the degeneracy between halo mass and environment 
for high mass haloes (see Figs, [l] and [2|. 

We have found strong and positive correlations in both 
the halo merger rate and the progenitor mass function with 
environmental densities. Figs. |4][7] present our main results, 
where haloes in the densest regions are seen to experience 
2 to 2.5 times higher merger rates than haloes in the voids. 
Such a density dependence can be approximated analyti- 
cally by multiplying our earlier fitting formula FM08 for the 
global merger rates (eqjsl by an additional (5-dependent fac- 
tor given by equation . This factor is a simple power-law 
in both the environmental density and halo mass, and it is 
redshift- independent. The mass dependence is quite weak, 
indicating that haloes with different masses but similar val- 
ues of 1 -I- St-fof experience similar merger rate amplifi- 
cations. This is intriguing in light of the fact, discussed in 
Section [3.2[ that these haloes actually reside in different en- 
vironments. 

The strong correlations of the halo merger rate and pro- 
genitor mass function with environment discussed in this 
paper have important implications for the analytic Press- 
Sch echter ([Press fc Schechter|1974|) a nd excursion set mod- 
els ( [Bond et al.||1991| [Lacey fc Cole||1993[ ). In this popular 
formalism, halo growth is modelled by the random walk tra- 
jectories of dark matter density perturbations smoothed at 
decreasing scales. Haloes are identified at scales at which 
these trajectories first cross some critical density threshold, 
and the Markovian nature of the model allows one to com- 
pute the distribution of these first crossings. This distribu- 
tion is then mapped onto the number-weighted conditional 
mass function (j>{M, z\Mo, zq) discussed in Section [4.2| and 
plays an important role in the Monte Carlo construction 
of mock merger tree catalogues (see Zhang et al.||200"8 and 
references therein). 

It is generally assumed that the conditional mass func- 
tion is independent of environment as the excursion set 
model is Markovian. The Markovian nature of the random 
walks, however, is not a prediction but rather an assump- 
tion resulting from the use of the fe-space tophat window 
function to smooth the density perturbations. There have 
been recent attempts to weaken this assumption or to in- 
troduce environmental dependence into other parts of the 
model (Zcntncr 2007, S andvik et al.|2007||Desjacquesp008l ), 
but these modifications thus far have not been able to repro- 
duce the basic statistical correlation between halo clustering 
and formation time found in simulation studies: older haloes 
are more clustered ([ Gottlober et al. ||2001 Sheth & Tormen 



2004 



2006 



Gao et al. 2005; Marker ct al. 2006; Wcchslcr et al 
Jing et al.. 2007,, Wang et al.,2007,,Gao fc White|2d07r 
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How do our environmental results for the merger rates 
tie in with these simulation and EPS studies? We have 
shown that the amplification of halo merger rates in denser 
regions persists at all redshifts (up to at least z = 2). If 
mergers were the dominant channel for halo growth, our 
results would imply that for haloes of a fixed mass today, 
those in denser regions should have formed more recently 
than those in void regions. Interestingly, this is exactly op- 
posite to the trend reported in many recent studies that have 
found older (i.e. earlier forming) haloes to be more clustered 
than younger haloes. As we will discuss in the next paper 
(Fakhouri & Ma 2008c), these two results are in fact not 



in confiict once the other important channel for halo mass 
growth - the "diffuse" accretion of non-halo material (ei- 
ther unresolved or stripped) - is taken into account. We will 
quantify the environmental dependence of this component 
and show that, when combined with the merger rate results 
presented in this paper, we recover the formation redshift 
dependence reported in prior simulation studies. 
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have unit area. The broadening of the 1 + <57_poF distribution 
with decreasing 2 is a natural consequence of gravitational insta- 
bility. This evolution is largely removed when ^{6j_-pO¥) is used 
as the variable, suggesting that i/'((57_fop) is a more self-similar 
measure of halo environment. 
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APPENDIX A: SELF-SIMILAR MASS AND 
ENVIRONMENT VARIABLES 

We have chosen to use the intuitive mass and environment 



variables M and Jt-fof in the fitting form (eq. Ill for the 
environmental dependence of the merger rate for ^ z ^ 
2. Here we investigate if other choices of mass and density 
variables may improve the fit. This is motivated by the well 
known property that the unconditional halo mass function 
is (almost) redshift-independent when the variable 



v{M,z) 



a{M)D{z) 



(Al) 



is used to characterise mass (see, for example, Jenkins et al. 
20011; whereas when M is used as the mass variable, the 



halo mass function evolves significantly with redshift. Here 
Sc = 1.686 is the critical overdensity for spherical collapse, 
a{M) is the variance of the linear density perturbations eval- 
uated at a scale corresponding to the halo mass M, and D{z) 
is the linear growth function. On the other hand, as we dis- 
cussed in the paper, the merger rate is more closely related 
to the conditional mass function than the unconditional mass 
function, and the redshift dependence of the former cannot 
be scaled out simply by using i/. Nonetheless, one can ask 
whether i/ is the more appropriate variable for capturing 
the mass dependence of the environmental dependence of 

dNmcigc/dz. 

A similar question can be raised about (Jy-pop. The 
overdensity S^-fof grows as a result of gravitational insta- 
bility, leading to broader distributions of 57_pof towards 
lower redshifts as shown in Fig. We can scale out the 
growth of S7-F0F in the linear regime by replacing ^t^fof (z) 
by 57-fof{z)/D{z). 

To incorporate the effects of nonlinear growth, we use 
the approach of [Goldberg & Vogeley ( 2004 1 for underdensi- 



ties and Peebles 



(19841 and Eke et al. (19961 for overdensl 



ties. These authors assume that the under/overdense regions 
are spherically symmetric and apply Birkhoff's theorem, 
treating these regions as self-contained universes embedded 
within an expanding universe. The model cosmological pa- 
rameters for these embedded cosmologies are computed from 
the density of the region under consideration, and the re- 
sulting Friedmann equation is solved to relate the densities 
at some redshift, 57^fof(z), to densities, i/j(57_pof(z)), at 
z = 0. When 57_fof ^ 1, this procedure is in agreement 
with the linear relation ^^{St-fof) ~ St-fof / D{z) 

A further complication for S^-fof is that the mass 
of the central object has been removed from 87. We have 
tested swapping the order of operation by first nonlinearly 
propagating 87 — > ip{S7) then subtracting the mass of the 
central object. We find that the resulting distributions of 
■!/;((57)_pop are only slightly modified from the distributions 

of ^(^7- FOf) . 

Fig. |A1| compares the distribution of the scaled 1 -f 
■!/;((S7_pof) (top panel) with that of the original 1 -I- 57_pof 
(bottom panel) for all haloes with M > 1.2 x lO^^M© at 
z = 0, 0.51, 1.08 and 2.07 (black to light grey lines). The up- 
per panel shows far less broadening with decreasing z than 
in the lower panel, indicating that -(/'(^t-fof) does remove 

^ We have also compared distributions across fixed v bins and 
found equivalent changes in the distribution of <57_pop. 



© 0000 RAS, MNRAS 000, [T]-[T7l 



,1.1 
i 1 

"0.9 



,1.1 

I 1 
"0.9 



, 1.1 

i 1 

"0.9 



L 1 ! 

-D.51 y- 


[M, 1 + dr-FOp] 


II III 1 1 1 MM 

-1.08 

...0. 


'l.08 ~ 

Q,5i 










1 1 1 Mill — 1 1 1 iiirt — ' 


i^im|^"'M''!^iiV-n 




-2.07 ' ' - 
...0. 




-2.07 ' ' - 
Q,51 

" 


-b.b7 - 
1,Q8 



Environmental Dependence of Halo Merger Rates 17 

[M, 1 + V(d>-FOF) 



,1.1 

i 1 



0.2 



1 10 0.2 1 10 0.2 

1 + '>~-mi- 1 + <57-FOF 



1 10 

1 + l57-FOF 




0.2 1 

1 + l/l(<57-FOF 



10 0.2 1 10 0.2 1 

1 1 + ^{Sy-FOF) 1 + *(<57-fof) 




0.2 1 10 0.2 

1 + 57-FOF 



10 0.2 1 10 

1 + 57-FOF 



-[M,S,z]/'- 



- [M, z] , computed 



1 + <57-FOF 

Figure A2. Ratios of the environmental dependence of the halo merger rates, R{M,S,z) — — — — 
at different redshifts: R{z2)/R{zi). Two mass and two environment variables are used, in clockwise order from the top-left matrix of 
subplots: (M, <57_pop),(M, 1 -|- i/'(<57-FOF)),('^(jW', z), 1 -|- '!/i(<57_pof))i (''(A^', z), (57_fof)- Each matrix presents ratios of R{z2) / R{zi) 
plotted as a function of environment. Each line represents a different mass bin with low mass in black and high mass in blue. The shaded 
regions represent Poisson errors. The redshifts used to compute R{z2) / R{zi) are noted in the top left corner of each subplot. The top-left 
matrix (M, <57_fof)i the variables used throughout this paper, shows the least systematic dependence on mass and environment. 



much of the redshift evolution of 57_fof- The mapping is 
clearly imperfect: The distributions at z = 1.08 and z — 2.07 
have long positive density tails that are not present at z = 0. 
This is not surprising since the simple spherical approxima- 
tion used for evolving the density cannot account for all 
non-linear effects such as the mergers of overdense and un- 
derdense regions. 

To test if 7/;((57_fof) and u are more appropriate vari- 
ables to use in the fitting formula than 5r-FOF and M, we 
show in Fig. A2 the ratio of dNmergc/dz[5] to the global 



^jdz for a number of log-spaced mass bins 
at 2 = 0,0.51,1.08, and 2.07. For brevity, let us refer to 
this ratio as H(M, 5, z). We can compute H using either 
M or i'(Af , z) as the mass variable, and either (57_fof or 
1-|-V'(57-fof) as the environmental variable. The upper left 
set of plots in Fig. |A2| uses M and 1 (57_fof, which are 
the variables used throughout this paper; the upper right set 
uses M and 1 -I- ■!/;(57_fof); the lower left set uses v(M,z') 
and 1 + 57_fof; and the lower right set uses v(M,z) and 
1 -I- V('57-fof)- 

Within each set of figures we plot a matrix of ratios 
of H computed at different redshifts. The redshifts used to 
compute the ratios are noted in the upper left corner of 
each subplot. For example, the upper subplot is the ratio 



H(M,<5,0.51)/H(Af,5,0) and is labelled "0.51/0". Each sub- 
plot contains five mass bins. Only points containing more 
than 40 haloes are plotted to minimise noise (this results in 
some mass bins being dropped). The variables that success- 
fully capture the redshift evolution in the merger rate will 
show very little variation in H(M, 5, z) with redshift and give 
ratios H(M, 5, zi)/H(M, 5, 22) close to unity. Interestingly, 
the [A/, (57_fof] pair used throughout this paper does the 
best job. The upper left matrix in Fig. |A2| shows ratios of 
H that tend to cluster around 1 and show few systematic 
trends with mass and environment. 

Using v{M, z) instead of M (lower panels) introduces 
a strong S dependence in the ratios of H: the 5-slope of 
dNmeige/dz flattens with increasing redshift. Similarly, using 
1 -I- ^{St-fof) instead of S^-fof (right panels) introduces 
a strong mass dependence when M is used as mass variable 
and does not improve on the 5 dependence introduced by 
u{M,z). 

Thus, Af and 1 + 57-fof appear to be the optimal vari- 
ables for capturing the environmental dependence of dark 
matter halo merger rates for ^ z ^ 2. 
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